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ABSTRACT 



We analyze Hubble Space Telescope surface-brightness profiles of 61 elliptical galaxies 
and spiral bulges (hereafter "hot" galaxies). The profiles are parameterized by break radius 
rb and break surface brightness These are combined with central velocity dispersions, 
total luminosities, rotation velocities, and isophote shapes to explore correlations among 
central and global properties. Luminous hot galaxies {My < —22) have cuspy cores with 
steep outer power-law profiles that break at r ~ v\, to shallow inner profiles / oc r -7 
with 7 < 0.3. Break radii and core luminosities for these objects are approximately 
proportional to effective radii and total luminosities. Scaling relations are presented for 
several core parameters as a function of total luminosity. Cores follow a fundamental plane 
that parallels the global fundamental plane for hot galaxies but is 30% thicker. Some of 
this extra thickness may be due to the effect of massive black holes (BHs) on central 
velocity dispersions. Faint hot galaxies (My > —20.5) show steep, largely featureless 
power-law profiles that lack cores. Measured values of r& and h for these galaxies are 
limits only. At a limiting radius of 10 pc, the centers of power-law galaxies are up to 
1000 times denser in mass and luminosity than the cores of large galaxies. At intermediate 
magnitudes (—22 < My < —20.5), core and power-law galaxies coexist, and there is a range 
in rb at a given luminosity of at least two orders of magnitude. Here, central properties 
correlate strongly with global rotation and shape: core galaxies tend to be boxy and slowly 
rotating, whereas power-law galaxies tend to be disky and rapidly rotating. A search for 
inner disks was conducted to test a claim in the literature, based on a smaller sample, that 
power laws originate from edge-on stellar disks. We find only limited evidence for such 
disks and believe that the difference between core and power-law profiles reflects a real 
difference in the spatial distribution of the luminous spheroidal component of the galaxy. 
The dense power-law centers of disky, rotating galaxies are consistent with their formation 
in gas-rich mergers. The parallel proposition, that cores are the by-products of gas-free 
stellar mergers, is less compelling for at least two reasons: (1) dissipationless hierarchical 
clustering does not appear to produce core profiles like those seen; (2) core galaxies accrete 
small, dense, gas-free galaxies at a rate sufficient to fill in their low-density cores if the 
satellites survived and sank to the center (whether the satellites survive is still an open 
question). An alternative model for core formation involves the orbital decay of massive 
black holes (BHs) that are accreted in mergers: the decaying BHs may heat and eject stars 
from the center, eroding a power law if any exists and scouring out a core. An average 
BH mass per spheroid of 0.002 times the stellar mass yields cores in fair agreement with 
observed cores and is consistent with the energetics of AGNs and the kinematic detection 
of BHs in nearby galaxies. An unresolved issue is why power-law galaxies also do not have 
cores if this process operates in all hot galaxies. 
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I. INTRODUCTION 



The Hubble Space Telescope (HST) allows us to study the centers of nearby galaxies 
with a resolution of a few parsecs. The centers of galaxies are interesting for several reasons: 
(1) some galaxy centers harbor AGNs and QSOs; (2) many or most galaxy centers may 
contain massive black holes that are the remnants of dead QSOs; (3) dynamical processes 
such as relaxation are more rapid near galaxy centers than elsewhere in the galaxy; thus 
interesting dynamical phenomena are likely to occur first near the center; (4) galaxy centers 
are to galactic astronomy as middens are to archaeologists: centers are the bottoms of 
potential wells and debris such as gas and dense stellar systems settle there, providing a 
record of the past history of the galaxy. 

The systematic properties of the centers of ellipticals and spiral bulges (hereafter 
"hot galaxies") were described by Lauer (1983, 1985a) and Kormendy (1982a, 1984, 1985, 
1987a,b). They detected inner regions in many galaxies where the slope of the surface- 
brightness profile flattens out, which they termed cores. They measured the size and 
surface brightness of these cores and demonstrated central parameter relations that linked 
core properties with one another and with global properties such as luminosity and effective 
radius. In particular, cores in brighter galaxies were larger and of lower density. The most 
recent version of the central parameter relations using ground-based data was presented 
by Kormendy and McClure (1993). A major goal of this paper is to revisit the central 
parameter relations using new iJSTdata on 61 galaxies. We shall show that HST broadly 
supports the ground-based scaling relations but elaborates upon them in important ways. 

Historically, the existence of cores in hot galaxies has been accepted as "normal" - 
probably because familiar dynamical models for galaxies such as the isothermal sphere and 
King models possess cores. In the absence of a central compact mass, it is plausible that 
all physical variables should vary smoothly near the origin and hence be expandable in 
a Taylor series with only even powers of r. In particular, the surface brightness may be 
written 

I(r)=I + hr 2 + O(r 4 ), (1) 

where r is projected radius. Using the conventional definition of core radius, a galaxy 
satisfying Eq. (1) would exhibit a core of radius r c such that I(r c ) = |/(0). Tremaine 
(1997) suggests the term "analytic core" for systems with cores in which all physical 
variables vary smoothly. The cores of King models and the isothermal sphere are thus 
analytic, while the i? 1 / 4 law is not. 

HST observations show that real cores are not analytic. In analytic cores, the surface 
brightness flattens at small radii as <i log // <i log r oc r 2 — note that this is stronger than 
the usual condition for a flat profile, d log I/d log r ^ — whereas real cores show shallow 
power-law cusps into the resolution limit (Crane et al. 1993; Kormendy et al. 1994; Jaffe 
et al. 1994; Lauer et al. 1995, hereafter Paper I; Kormendy et al. 1996a). Fits in Byun 
et al. (1996, Paper II) yield projected slopes 7 = — <i log i"/<i log r in the range 0.05-0.3 for 
surface brightness, while non-parametric inversions for space density show even steeper 
slopes, from 0.2 to 1.1 (Gebhardt et al. 1996, Paper III; and Kormendy et al. 1996a). 
Thus, real cores have divergent rather than constant densities as r — > 0. 
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So far, cores have been found only in luminous ellipticals. The division between core 
and non-core galaxies is fairly sharp. Surface-brightness profiles either flatten out to form 
cores or continue to rise steeply into the resolution limit — few galaxies are in between 
(Kormendy et al. 1994; Jaffe et al. 1994; Paper I; Kormendy et al. 1996a). Statistical 
analysis of non-parametrically derived space density profiles indicates the existence of two 
groups (core and non-core) at the 90% confidence level (Paper III). 

Paper I introduced the term power laws to describe the steeply rising, featureless 
profiles that lack cores. 5 It is possible that the power-law category as we have drawn 
it may be oversimplified: At present the category contains a number of low-luminosity 
galaxies whose upper limits on core size are larger than those predicted by extrapolation 
of the core-luminosity relationship defined by brighter galaxies. In other words, some of 
the low-luminosity power-law galaxies may really be part of a core sequence extending 
to lower luminosity. Recent WFPC2 images in fact show tiny cores in a few power-law 
galaxies (Lauer et al. 1997). Nevertheless, the upper limits on core size for brighter power- 
law galaxies are already well below the core sequence for galaxies of similar luminosity, and 
thus clearly differentiate them. Future results may compel some revision of the power-law 
category, but the present simple core/power-law division is a useful working hypothesis. 

Lauer (1985a) emphasized that the central properties of hot galaxies do not correlate 
perfectly with total luminosity and sought an explanation in terms of a second parame- 
ter. The present data suggest that this second parameter is related to global rotation and 
isophote shape. So far, cores have been found only in luminous, slowly rotating ellipticals 
with boxy isophotes 6 , while power laws are found in faint, rapidly rotating galaxies with 
disky isophotes. A link between central profile type and global shape/rotation was sug- 
gested by Nieto et al. (1991a) based on ground-based images, and further evidence was 
presented by Jaffe et al. (1994) and Ferrarese et al. (1994) based on HST images of 14 Virgo 
galaxies. The present database is considerably larger and permits a critical examination 
of this link and its relation to hot galaxy formation. Our point of view differs importantly 
from that of Jaffe et al., who ascribe many of the differences between the two profile types 
to inclination effects connected with a small inner disk seen either face-on or edge-on. 
In contrast, we — like Nieto et al. - believe that the spheroidal light distributions are 
intrinsically different in the two types and would look the same from any viewing angle. 
These differences in viewpoint are discussed in Section 5 and Appendix A. 

The results we have described raise several theoretical issues: why are there two types 
of profile and how did each type form? Why do the two types have different global rotation 
and shape? Why are cores non-analytic? And what do central profiles tell us about hot 
galaxy formation and evolution? 

The second, more speculative, part of this paper addresses these issues. We suggest in 
Section 7 that the power-law profiles of disky galaxies indicate dissipation and are therefore 
consistent with formation in gas-rich mergers. The parallel suggestion — that the cores of 

5 Cores and power laws were also identified by Jaffe et al. (1994), who called them Type 
I and Type II. 

6 Isophote shape in elliptical galaxies is explained and defined by Bender & Mollenhoff 
(1987). A recent discussion is given by Kormendy & Bender (1996). 



- 5 - 



boxy galaxies are the by-products of purely stellar, gas-poor mergers — is more problematic. 
For example, luminous core galaxies are expected to accrete small dense satellites. The 
rate of such accretions appears sufficient to gradually fill in all low-density cores if such 
satellites survived and sank to the center. An unresolved issue is whether the satellites do 
survive, and thus whether some other process is needed to defend low-density cores against 
in-fill. 

Even if the data do not firmly require such a mechanism, there is strong and growing 
evidence for a widespread population of massive central black holes (BHs) in hot galaxies 
(Kormendy & Richstone 1995). The presence of these objects must be taken into account 
in standard merger-based models for forming hot galaxies (Section 8). The BHs associated 
with the merging galaxies form binaries whose orbits then decay. The orbital decay heats 
the surrounding stars, erodes a power law if one exists, and scours out a core. Accreted 
satellites will also tend to be ripped apart, thus preventing core in-fill. BHs with plausible 
masses (as estimated in Appendix B) seem able to produce cores of roughly the right 
size and scaling versus galaxy luminosity. In this way, the presence of central BHs might 
"rescue" the dissipationless, gas-poor model for cores and boxy galaxies. However, models 
of core formation based purely on massive BHs leave several questions open, notably how 
power-law profiles escape similar disruption. 

Whether or not these speculations about galaxy formation are correct, the updated 
relations between central and global galaxy parameters that are presented in this paper 
appear to provide important new constraints on hot galaxy formation. 

2. CENTRAL PROFILE TYPES 

Major collections of HST central profiles include Crane et al. (1993), Jaffe et al. (1994), 
Forbes et al. (1995), and Paper I. An assortment of representative surface- brightness 
profiles of 55 ellipticals and spiral bulges is given in Fig. 1. The following summary is 
based on the data and discussion in Paper I. 

We distinguish two types of hot galaxy: 

(1) Core galaxies have "broken" power-law surface-brightness profiles that change slope 
significantly at a "break radius" r&. To identify a galaxy as having a core, we re- 
quire that the absolute value of the inner logarithmic slope, 7 = — <i log i"/<i log r, be 
shallower than 0.3. Nearly all core galaxies appear to have 7 > 0, which indicates 
a cusp in the central surface brightness and an even stronger cusp in the luminosity 
density. Paper III concluded that, even with errors taken into account, only 2 out of 
15 known core galaxies could admit an analytic core (7 = 0). Core galaxies as a class 
are luminous objects with My < — 20.5 (H = 80 km s -1 Mpc -1 ). They range from 
brightest cluster galaxies down to the intermediate-mass field elliptical NGC 3379. 

(2) Power-law galaxies show fairly steep surface-brightness profiles with no significant 
break within 10" (at Virgo). Their average surface- brightness slope is 7 ~ 0.8 ± 0.3 
at the smallest resolvable radius. Power-law galaxies are generally fainter than core 
galaxies (My > —22), but their luminosity densities at 10 pc are 10-1000 times higher 
than those of cores (Paper I). Profile shapes within .1 are generally not known, 
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though recent WFPC2 images suggest small cores inside some power laws. Power-law 
galaxies include M 32 (NGC 221), small Virgo ellipticals, and bulges of disk galaxies. 
Both profile types are well fit by the following equation (the "Nuker" law, Papers I 



The asymptotic logarithmic slope inside r\, is —7, the asymptotic outer slope is — /?, and 
the parameter a parameterizes the sharpness of the break. The break radius r& is the 
point of maximum curvature in log-log coordinates. "Break surface brightness," Ij,, is the 
surface brightness at r&. Equation (2) is intended to fit only over radii accessible to the 
HST Planetary Camera, i.e., <10". For typical fitted values of /?, there must be a further 
turndown in the profile at larger radii for the total luminosity to be finite. 

Nuclei are identified when excess light above the prediction of Eq. (2) is visible within 
the inner few tenths of an arcsec. Nuclei with varying degrees of prominence are illustrated 
in Paper I (Fig. 14). Objects with prominent nuclei are always systems of low luminosity 
and are probably nucleated dSph or dE galaxies. Nuclei are assumed to be star clusters 
(or possibly unresolved tiny stellar disks), but direct spectral confirmation is often lacking. 
A stellar nucleus in NGC 3115 has been resolved in recent WFPC2 images (Kormendy 
et al. 1996b). Non-thermal central point sources exist in four galaxies in our sample: M 87 
(NGC 4486), NGC 6166, Abell 2052 (Paper I) and NGC 4594 (Kormendy et al. 1996c). 
We call these AGNs to distinguish them from nuclei. So far, no nuclei (as opposed to 
AGNs) have been found within cores (Kormendy & Djorgovski 1989; Paper I). 

Resolution plays an important role in classifying profiles and estimating central prop- 
erties. This is illustrated in Fig. 2, which shows M 31 (NGC 224) and M 32 (NGC 221) as 
seen at their actual distances and as they would be seen 24 times further away just beyond 
Virgo (for future reference, we call these artificially positioned galaxies M 31-in- Virgo and 
M 32-in-Virgo) . Up close, M 31 shows a two-component profile that is clearly divided into 
a bulge and a nucleus, the latter showing a small core. The entire profile shows too much 
substructure to fit comfortably into either the core or power-law category. In contrast, 
M 31-in- Virgo shows only a trace of a nucleus, and its profile and degree of nucleation 
are similar to those of several other galaxies that we have classed as power laws (see Fig. 
14 in Paper I for a collection of power laws with varying degrees of nucleation). M 31 
implies that many power-law galaxies, particularly those with hints of nuclei, may contain 
significant substructure, including nuclei and tiny cores. 

M 32 is similarly ambiguous. Seen up close, M 32's profile in Fig. 2 breaks from a 
power law near .5, curving gently downward into the resolution limit. M 32-in-Virgo 
shows a nearly perfect power law with only a small bend at the equivalent nearby radius of 
70". Thus M 32 does not fit Eq. (2) very well either, because of multiple breaks that yield 
different values of rj, depending on what portion of the profile is fitted. M 32 shows that 
values of r& in power-law galaxies are not robust and that similar breaks at small radii 
could exist in other distant power-law galaxies, even those that apparently show clean 
power laws at the present resolution. 
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Because the fitted values of r& in power-law galaxies are less robust than those for 
core profiles, which reflect real features, we regard them as less fundamental. As explained 
below, we treat the fitted values of r& differently in analyzing the two types of galaxy. 

3. GALAXY SAMPLE AND DATABASE 

The database used in this paper is contained in Tables 1, 2, and 3. A brief overview is 
given here, and additional details are provided in the table notes. The heart of the sample 
consists of 42 normal ellipticals and bulges taken from Paper I (NGC 4150, NGC 4826, and 
NGC 5322 were excluded due to strong nuclear dust). To these were added images of 14 
E's and bulges from the WFPC1 GTO programs (some unpublished). Five more normal 
E's, mostly Virgo galaxies from Jaffe et a/., were located in the HST public archive as of 
June 1993, for a total of 61 galaxies. The original GO/GTO program and references to 
published HST profiles are listed in Table 1. All images were taken using the Planetary 
Camera in Cycles 1 and 2 and consequently suffer spherical aberration. They were observed 
through filter F555W, which approximates the V band, and usually have a peak signal of 
> 10 4 photons in the central pixel. All images were processed as described in Paper I and 
deconvolved with the same Lucy-Richardson procedure used there. 

Power-law galaxies with identified nuclei are divided into two types: "moderately" and 
"severely" nucleated, indicated in Table 1 by "+" and "++"• M 31-in- Virgo is adopted 
as the dividing line between the two types (cf. Fig. 2 here and Fig. 14 of Paper I). In 
severely nucleated galaxies and in galaxies with AGNs, fits to the nuker law ignore the 
innermost pixels affected by the nuclear light. 

Table 1 presents observed quantities such as Hubble type, distance, magnitude, color, 
and nuker-law parameters from Paper II. A few galaxies not treated in Paper II have been 
similarly fit and the results are given here. M 31 and M 32 appear twice, as seen nearby 
and near Virgo (labeled with a "V"). For core-type profiles, we accept the nuker-law fits as 
given for 9b and fib- 7 For power-law galaxies, no core is resolved, and we use the separate 
upper limits on core size and surface brightness provided by Paper I. These limits (for 
power laws only) are called 9 l h im and (J^ m in Table 1. For a few power-law galaxies not 
contained in Paper I, these limits were obtained from a visual estimate of the steepness of 
the innermost part of the profile. 

The distance to each galaxy (in km s _1 ) has been estimated using a variety of methods 
as summarized in the notes to Table 1, and the adopted value and its conversion to Mpc 
(based on Hq = 80 km s _1 Mpc _1 ) are given there. These distances are used to convert 
the apparent quantities in Table 1 to absolute quantities in Table 2. B-band magnitudes 
are converted to the V band to be consistent with the HST profiles. Data taken from 
the literature include central velocity dispersion, o"o, an inner velocity dispersion gradient 
defined as R a = ao/a(10"), dimensionless rotation parameter (t>/cr)*, isophote shape a^/a, 
global (effective) radius r e , and global surface brightness /U e , defined as the mean surface 
brightness within r e . Details and sources are given in the notes. 

7 The quantity 9b is the break radius in arcsec, while fib is break surface brightness 
expressed as V mag arcsec -2 . 



- 8 - 



Table 3 presents several derived quantities based on spherical, isotropic dynamical 
models fitted to the nuker-law light profile. The mass-to-light ratio of each model has 
been determined by normalizing to o~q from Table 1, assuming constant M/L with radius 
and equating ao to the light-weighted rms line-of-sight dispersion in a centered 2" by 2" 
aperture (corrections for 1" FWHM seeing are at most a few percent and are not included). 
Mass-related quantities are blank if gq is not available. Quantities tabulated at .1 include 
the luminosity density, peak Maxwellian phase-space density, two-body relaxation time, 
and predicted projected velocity dispersion. Total luminosity and mass within a sphere 
of the same radius are also given. Comparison to the non-parametric densities in Paper 
III indicates that nuker-law fitted luminosity densities are 10% too low on average but 
otherwise show little scatter for non- and moderately nucleated galaxies (severe nuclei 
were ignored in fitting nuker laws, and as a result nuker-law densities in these galaxies 
are about a factor of 2 lower than the non-parametric inversions). Several quantities are 
repeated for r = 10 pc, but for many galaxies this is well inside the resolution limit of .1 
and requires an inward extrapolation of the nuker-law fit. 

An impression of the division into core and power-law galaxies is provided by Fig. 3, 
which plots inner power-law slope 7 versus observed break radius 6b (or 6 l b im for power-law 
galaxies) in arcsec. Profiles with 6b > .16 (log 6b > —0.8) are reasonably well resolved 
by HST. They divide into two groups, those with 7 < 0.25 (cores) and those with 7 > 0.5 
(power laws) — the valley in between is empty. This is the division that motivated the 
two profile types in Paper I, later analyzed statistically in Paper III. 

The rectangular box in Fig. 3 encloses galaxies that we are fairly sure contain real 
cores. Galaxies above the box are definitely power laws at current resolution. Galaxies to 
the left of the box are classed as power laws although some contain a hint of an incipient 
core. The effect of limited resolution is apparent for M 31 and M 32; both galaxies are 
plotted twice, as seen nearby and at Virgo. The plotted positions differ appreciably, 
reflecting features of their inner profiles that cannot be probed in more distant galaxies. 

Galaxies within the box in Fig. 3 comprise the "Core" sample used in the following 
section. All others are classed as power laws. 

4. CENTRAL PARAMETER RELATIONS 

The data in Tables 1 and 2 are used to plot new central parameter diagrams like 
those of Lauer (1983, 1985a) and Kormendy (1985, 1987a,b). We begin with plots versus 
absolute magnitude in Figs. 4a,b,c,d. The symbols have the following meanings: 

(1) Core galaxies are plotted with filled circles (•) using values of r& and from Table 2. 

(2) Power laws are plotted with open circles (o) using the limits r l h %m and fi l b im from Table 
2. 

(3) M 31 and M 32 are plotted twice, as seen at their actual distance (asterisks) and in 
Virgo (end of vector). The length and direction of these vectors illustrate the possible 
effect of changing resolution on other power-law galaxies. Their direction is opposite 
to the limit flags that are attached to all power-law galaxies. 
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(4) Special objects: The SO galaxy NGC 524 is the only core profile that is found within 
a bulge (all others are in ellipticals). NGC 524 is roughly face-on and shows flocculent 
dusty disk arms (Paper I) and a blue center (Kormendy, private communication); it 
is plotted with a small square. Fornax A (NGC 1316) is a probable recent merger 
remnant (Schweizer 1980) with a peculiar morphology (RC3). It has an abnormally 
small core for a galaxy of its luminosity (Kormendy 1987b). NGC 4486B shows a 
double nucleus like M 31 's in WFPC2 images (Lauer et al. 1996) but continues to 
have a clearly defined core. 

The new plots show the same broad trends versus galaxy luminosity that were seen in 
ground-based data (Kormendy & McClure 1993). Core galaxies are luminous objects that 
extend down to My = —20.5. All normal ellipticals brighter than My = —22 show cores, 
with cores of brighter galaxies being larger and lower in surface brightness and density. 
The new central parameters of core galaxies correlate well with previous values measured 
from the ground (Kormendy et al. 1994). The parameter relations for core galaxies are 
fairly narrow; for example, the rms scatter in r& versus My about the best-fitting line is 
only 0.25 dex (Fornax A omitted). 

Ferrarese et al. (1994) have questioned whether the trends in core properties versus 
absolute magnitude are an artifact created by adding brightest cluster galaxies (BCGs) 
to smaller, trendless galaxies. They argue that, aside from M 87, all cores in their Virgo 
sample are of similar size, and trends appear only when M 87 is added. Although M 87 
does not strictly qualify as a BCG (that distinction in Virgo is held by NGC 4472), it does 
share certain properties with BCGs such as high luminosity and central location within a 
sub cluster. 

From our larger sample, it seems clear that trends in core properties versus My are 
real and are not an artifact of adding BCGs. The present sample could be truncated 
at My = —22.2 to eliminate all BCGs (including those in small groups as well as Abell 
clusters), yet trends among the 11 remaining core galaxies between My = —20.5 and —22 
would still be present. In all plots, core properties of BCG galaxies appear to be a normal 
extension of the cores in smaller core ellipticals. 

Power-law galaxies in Fig. 4 are low-to-intermediate luminosity systems that extend in 
luminosity up to My = —22. They overlap with core galaxies at intermediate magnitudes 
in the range —20.5 > My > —22. Despite an increase in angular resolution by a factor 
of 10 with HST, we have generally failed to find cores in these objects, and thus their 
distribution in Fig. 4a is rather flat, reflecting the constant HST resolution limit of ~ 0.1 
arcsec. For systems fainter than My —19, this limit is uninteresting since it equals or 
exceeds predictions based on extrapolation from core galaxies. However, at intermediate 
magnitudes in the range My = —20.5 to —22, power-law and core galaxies coexist, and 
it is clear that the scatter in break radius is real and large. Core/power-law pairs that 
illustrate extremes of r\, at fixed luminosity include NGC 3379 and NGC 1023, whose break 
radii differ by more than a factor of 40 while their absolute magnitudes differ by less than 
0.5 mag, and NGC 4168 and NGC 4594, for which the ratio of break radii is over 100 
even though their absolute magnitudes are almost identical. This is not a resolution effect 
wherein cores are detected in nearby galaxies but not in distant ones. Figure 5 plots break 
radius versus distance and shows that most of the sample, containing both small and large 
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cores, resides in a narrow range of distance near that of Virgo. More distant galaxies are 
actually more likely to show cores because their cores are intrinsically larger. 

The large scatter in break radii near My = —20.5 to —22 might at first sight be taken 
as a manifestation of the two-dimensional, planar distribution of the global structural 
parameters of hot galaxies, i.e., the fundamental plane (Dressier et al. 1987; Djorgovski 
& Davis 1987; Faber et al. 1987). Two-coordinate projections of this two-dimensional 
distribution commonly exhibit scatter depending on whether they show the plane edge-on 
or face-on. The basic coordinates for the global plane (see Section 6) are r e , \i e , and do, 
from which Ly can be derived as Ly = 27Tfx e r^. A plot of radius versus magnitude is thus 
a projection of the fundamental plane, and scatter might be expected in r& versus My that 
is comparable to that seen in r e versus My, provided and r e are well correlated. 

This hypothesis is tested by substituting r e for r b in Fig. 4d. The scatter there proves 
to be small, demonstrating that the combination of radius versus L shows the global plane 
rather close to edge-on. The much larger scatter of Fig. 4a therefore suggests a real 
decoupling of central properties from global ones, as emphasized by Lauer (1985a). In 
Section 5 we examine this scatter in more detail and show that it correlates with global 
rotation and isophote shape, in the sense that power-law galaxies (which have small r&) 
are disky and rotate rapidly, while cores (which have large r&) are boxy and rotate slowly. 

Bulges are distributed in Fig. 4 like ellipticals of small-to-intermediate size. None 
(except for M 31-nearby) shows a core. The resemblance of bulges to small and intermediate 
ellipticals is not surprising since the two classes of galaxy share several traits, including 
similar global size, high rotation, flattening by rotation rather than anisotropy, and disky 
subsystems (Bender, Burstein & Faber 1992). 

Before drawing further conclusions from Fig. 4, we consider whether the trends shown 
there are affected by the particular sample of galaxies chosen. The present sample is a 
mixture taken from different authors, but we have been careful to retain only objects that 
are morphologically normal and free of dust. Our own sample from Paper I (comprising 
42 out of the 61 total objects in this paper) was specifically chosen to probe the full 
range of parameters covered by the ground-based central parameter relations (Lauer 1985a; 
Kormendy & McClure 1993). We strove hard to sample the widest possible magnitude 
range and, at intermediate magnitudes, to sample galaxies with both large and small 
apparent cores. Thus, it is possible that the present sample somewhat exaggerates the 
total spread in break radii at middle magnitudes. 

Another point is that most objects studied here had previous ground-based data, and 
thus some prior clue as to core size. Since ground data typically agree well with HST data 
(especially for large galaxies, Kormendy et al. 1994), the present sample does not provide 
a truly fresh look at galaxy centers. A sample to do this with completely new galaxies 
has been observed in Cycle 5 and is now being analyzed. What the present sample does is 
fairly probe galaxies that had previously been examined from the ground. 

Are the claimed correlations robust for core galaxies specifically? Although Fornax 
A has been included in the diagrams for interest, it is strongly peculiar and its center is 
contaminated by dust (Shaya et al. 1996). It does not qualify for our sample of normal, 
massive E's, and its high residuals should not count against the correlations. Six more 
core ellipticals with ground-based data could also be added to bolster the HST data; these 
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new galaxies agree well with the trends here (Kormendy, unpublished). Thus we feel that 
evidence strongly favors the core correlations found in Fig. 4; however, full confirmation 
will require the completely independent sample of galaxies from Cycle 5 that we are now 
analyzing. 

Finally, the number of galaxies at intermediate magnitudes is still small. We will argue 
in the next section that the spread of properties in this magnitude range is correlated with 
global boxiness and diskiness, and will draw from this some significant conclusions about 
hot galaxy formation. Clearly, these conclusions will need to be checked by gathering a 
larger and more representative sample of galaxies at these magnitudes. Again, the Cycle 
5 sample was selected to do this. 

The referee has asked whether bright power-law galaxies might in fact all be SO's (or 
bulges). There are 7 power-law galaxies in the magnitude range —20.5 to —22: NGC 596 
(E+4:), NGC 1172 (E+2:), NGC 1700 (E4), NGC 3115 (SO"), NGC 4594 (Sa), NGC 4621 
(E5), and NGC 4697 (E6) (Hubble types from the RC3). Two of these are actual disk galax- 
ies (NGC 3115, NGC 4594), two more are highly flattened E's (NGC 4621, NGC 4697), 
and two more have SO-like outer envelopes (NGC 596, NGC 1172). That leaves only 
NGC 1700, which is also fairly flattened. We will show in the next section that there 
is a good correlation between power-law galaxies and galaxies that are rapidly rotating 
with disky isophotes. Thus it is possible that all bright E-type power-laws are in fact SO's 
masquerading as ellipticals. 

It is interesting to speculate where the bright power-law galaxies will move within 
Fig. 4 as resolution improves. At present these galaxies fall below the core sequence by 
x3-10 in core size. However, new WFPC2 observations have increased this distance for a 
number of them (Lauer et al. 1997). It is thus possible that bright power-law galaxies may 
ultimately emerge as a separate population rather than simply tracing the lower edge of a 
large dispersion in core size at these luminosities. 

M 31 and M 32 again show the effect of changing spatial resolution. M 32-in-Virgo lies 
near the Virgo dwarf E's in all of its parameters, suggesting that the downward trend in 
surface brightness for faint galaxies in Fig. 4c is mainly an artifact of resolution. M 31-in- 
Virgo is indistinguishable from other bulges of similar magnitude. Its shift in the diagrams 
under distance change is not as large as that of M 32 because its profile is not as steep as 
M 32's between O'.'l and 2 "4 (Fig. 2). 

NGC 4486B is the lone core galaxy near My = —17.5 mag. Its angular size lies near 
the lower boundary of established cores in Fig. 3, but the presence of a core has been 
confirmed in WFPC2 images (Lauer et al. 1996). NGC 4486B's low luminosity, compact 
profile, high line-strength, and close proximity to M 87 suggest that it might be tidally 
stripped by its larger neighbor (Faber 1973). Its core parameters would be consistent if it 
once resembled the small-core galaxy NGC 3608 and then lost ~90% of its outer luminous 
envelope. Whether a core could actually survive such extensive stripping and whether a 
diffuse giant like M 87 could strip a high-density object like NGC 4486B are open questions 
(see Section 7.3.2). 

Figure 6 plots luminosity density, mass density, and phase-space density at a limiting 
radius of 0.1 (from Table 3) as a function of absolute magnitude. An impressive feature 
of Fig. 6 is the large range in density from small to large galaxies, almost 10 6 in all three 
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panels; half of this range is spanned by cores alone. Part of this spread is due to the fact 
that Fig. 6 mixes objects at different distances. However, densities at the fixed physical 
scale of 10 pc still show a range of 1000 (column 13, Table 3; see also Fig. 13, Paper I; 
note that densities in distant galaxies at 10 pc require inward extrapolation). The large 
range of densities near the centers of hot galaxies has been remarked on before (Kormendy 
1984, Lauer 1985a, Carlberg 1986), but HSTs higher resolution has pushed up densities 
in power-law galaxies by another two orders of magnitude. As discussed in Section 7, the 
large density contrast between small and large hot galaxies is an important clue to hot 
galaxy formation. 

We conclude this section by presenting general scaling laws for core properties versus 
galaxy luminosity and mass. We restrict attention to core galaxies because their param- 
eters are robust. It is well known (see Section 6) that cores are really a two-dimensional 
dynamical family (the fundamental plane). Nevertheless, it is often convenient to treat 
them as a one-parameter family depending on luminosity or mass. This is possible be- 
cause (1) the fundamental plane is only partially filled, and projections against any pair 
of coordinate axes have only limited extent; (2) typically, these projections are elongated 
and can be approximated by one-dimensional scaling relations; and (3) core mass and lu- 
minosity are both moderately well related to global mass and luminosity for core galaxies. 
This combination produces the tight correlations versus total luminosity in Figs. 4 and 6, 
which we now quantify. 

The following is a set of self-consistent scaling relations versus Ly and galaxy mass, 
M. Total luminosity has been related to mass by assuming (M/Ly) oc L 25 (Faber 
et al. 1987). The exponents in these relations are not least-square fits but have been 
derived by a process of trial-and-error adjustment to maintain consistency with standard 
structural formulae. The first three of these relations are independent fits to the data, 
while the rest are derived from the structural formulae. The relations involving M core 
assume that the core is in dynamical equilibrium. 
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5. SCATTER IN THE CENTRAL PARAMETER RELATIONS: 
CORRELATIONS WITH GLOBAL ROTATION AND ISOPHOTE SHAPE 



An emerging suspicion of the last decade is that there are actually two types of elliptical 
galaxies: luminous E's with boxy isophotes that rotate slowly, and small E's with disky 
isophotes that rotate rapidly (Bender 1988; Nieto 1988; Nieto & Bender 1989; Bender 
et al. 1989; Kormendy & Djorgovski 1989). A formal division of the Hubble sequence for 
ellipticals into two classes has been suggested based on these criteria (Kormendy & Bender 
1996). We shall refer to these two subtypes as boxy and disky respectively. 

The scatter seen in central parameters versus absolute magnitude in Figs. 4 and 6 
appears to correlate with boxy/disky subtype. This correlation is illustrated in Fig. 7, 
which replots Fig. 4a (r b versus My) but now with symbols indicating rotation, (u/cr)*, 
in Fig. 7a, and isophote shape, a^/a, in Fig. 7b. It is seen that power-law galaxies are 
mainly rapidly rotating and disky, while cores are slowly rotating and boxy or neutral. At 
intermediate magnitudes, My = —20.5 to —22, the presence of a core is a better predictor 
of boxiness or slow rotation than is absolute magnitude. The same correlation was also 
found by Nieto et al. (1991a) using ground-based data. In our sample, the correlation 
between disky galaxies and power laws appears strongest if the disky threshold is set at 
a^/a > 0.4. All other isophote types (boxy, neutral, and variable) are associated with 
cores. "Variable" galaxies are those with strongly varying values of a^/a versus radius. 
We use the term boxy to include the boxy, neutral, and variable types. 

Jaffe et al. (1994) and Ferrarese et al. (1994) also saw a link between power-laws and 
disks based on their HST Virgo sample. They went further to suggest that essentially all 
power laws have inner disks within 1" and that such disks seen edge-on are what produce 
power-law profiles (Jaffe et al. 1994). We agree that power laws tend to reside in galaxies 
with globally high rotation and global diskiness, but we do not see evidence for inner disks 
in all or even most power laws. Rather, we believe that the hot component is intrinsically 
different in cores and power laws and that the two classes therefore would look different 
from any viewing angle. 

This difference in interpretation is fundamental, since implications for galaxy forma- 
tion would be limited if profiles were a sensitive function of viewing angle. We believe 
that Jaffe et al. were influenced by the fact that most power-law objects in their sample 
happened to be highly flattened, edge-on, late- type E/SO's with a higher-than-average in- 
cidence of both inner and outer disks. Our power-law galaxies as a group are less flattened, 
less edge-on, less skewed to late Hubble types, and do not in general show inner or outer 
disks. A brief discussion of the edge-on disk model for power laws was presented in Paper 
I. A more extensive search for inner disks was made with the present larger sample, and 
the results are reported in Appendix A. 

To summarize, it appears likely that disky and boxy ellipticals have different kinds of 
central light profiles. Since global properties are implicated, whatever process established 
this connection was probably a major event in the life of the galaxy. We explore this link 
and its implications for hot galaxy formation in Section 7. 



-14- 



6. THE CORE FUNDAMENTAL PLANE 



So far we have focussed on the relationship between the centers of galaxies and their 
global properties. We turn next to relationships among the central properties alone. Plot- 
ting galaxies in central (log rb, fj,},, logo" )-space allows us to look for a fundamental plane 
(FP) analogous to the one found in global (logr e , jU e , log cr e )-space (Dressier et al. 1987; 
Djorgovski & Davis 1987; Faber et al. 1987). If cores are in dynamical equilibrium (highly 
likely), if they are supported by random motions (as indicated by the observations), if ve- 
locity anisotropy does not vary too much from galaxy to galaxy (unknown, but see below), 
and if core M/L is a well-behaved function of any two variables r&, or o"o (true of global 
M/L), then we expect cores to populate a thin surface in central (logr^, Hb, logo"o)-space 
(Faber et al. 1987). Lauer (1985a) demonstrated that cores of well-resolved galaxies indeed 
populate a two-dimensional manifold, and Faber et al. (1987), using Lauer's data, derived 
a preliminary core fundamental plane that was roughly parallel to the global FP. We revisit 
the tilt and thickness of this plane using the more accurate HST data. 

Figure 8 shows projections of (logr;,, /if,, loga )-space for the present sample. Cores 
are again filled circles, while power-law galaxies (open circles) have been plotted for com- 
pleteness' sake using their limiting values (warning: the direction of the limit flags is 
schematic using M 31 and M 32 as a guide). To seek a plane, we rotate about an axis and 
search for the thinnest distribution of points. Following Faber et al. (1987), we choose to 
rotate about the ao-axis because the resulting combination of lb and is nearly invariant 
to small resolution changes and small measurement errors. The best rotation (based on 
core galaxies only) is shown in Fig. 8d. Within the errors, the tilt of this FP is consistent 
with the global plane, a 135 ~ r e /g 84 , found by Faber et al. (1987); this is the orientation 
plotted. The residual rms scatter about this plane, expressed as an error in logr^, is 0.12 
dex (for cores only). This is 30% larger than the equivalent scatter about the global FP, 
which is 0.09 dex (Lynden-Bell et al. 1988). The larger scatter about the core plane may 
be related to the presence of central BHs, which elevate o"o in some galaxies. Five labeled 
objects with positive residuals in Fig. 8d are BH candidates from Kormendy & Richstone 
(1995) and Kormendy et al. (1996b, c). 

The locus of h versus rj, is also quite narrow for cores (Fig. 8c). This occurs because 
the individual brightness profiles of core galaxies are approximately tangent to a single 
line in the (logr?,, log/;,) plane, as can be seen in Fig. 1. This tight correlation can be 
used to construct a second distance indicator based on and lb alone. The scatter about 
the best-fitting line corresponds to 0.18 dex in log distance, which is two times worse than 
the global FP. However, the method does not require any measurement of a and thus may 
sometimes be useful. 

The existence of a core fundamental plane suggests that (1) cores are in dynamical 
equilibrium supported by random motions; (2) that r& and lb are meaningful dynamical pa- 
rameters describing the size and luminosity density of the core; (3) that velocity anisotropy 
does not vary greatly among core galaxies; (4) that the mass of any central BH does not 
strongly dominate the core potential in most galaxies; and (5) that core mass-to-light ratio 
varies smoothly over the fundamental plane. 
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The core fundamental plane is well defined even though (1) the profiles of core galaxies 
are not analytic; (2) the profiles of different core galaxies are not identical (a and 7 vary, 
Paper I); (3) the velocity dispersion anisotropy may vary from galaxy to galaxy; (4) some 
or all core galaxies may harbor massive BHs that distort both the photometric profile 
and do; and (5) cores represent only a tiny fraction of the total luminosity of the galaxy. 
Evidently whatever differences exist among core galaxies are not so large as to erase the 
appearance of a two-parameter family of self-gravitating cores that is fundamentally not 
too dissimilar from the 2-dimensional family of isothermal spheres. 

7. CENTRAL PARAMETER RELATIONS AND HOT GALAXY FORMATION 

The final sections of this paper discuss the central parameter relations in the context of 
galaxy formation. We assume throughout that hot galaxies form via hierarchical clustering 
and merging (hereafter HCM; see e.g., Toomre 1977, White & Rees 1978, and Blumenthal 
et al. 1984). Descriptions of HCM as applied to hot galaxies may be found in Schweizer 
(1986), Kormendy & Sanders (1992), de Zeeuw & Franx (1991), Barnes & Hernquist (1992), 
and Bender et al. (1992). An important challenge to HCM is the formation of boxy and 
disky galaxies, including the association found here with central cores and power laws. A 
novel element that needs to be considered is the presence of massive BHs in the centers 
of many or most hot galaxies, which could considerably alter the predictions of standard 
HCM with no BHs. 

We will suggest that the high-density cusps of power-law galaxies are broadly consis- 
tent with the growing body of evidence that points to the importance of gaseous dissipation 
in their formation. More of a puzzle are the low-density cores of massive hot galaxies - 
they seem hard to form in the first place and hard to maintain once formed. The problem 
of cores leads us to consider an alternative method for forming them based on central BHs. 
This model is the subject of Section 8. 

7.1 Previous Work on Hot Centers 
Many authors have discussed the centers of hot galaxies in the context of galaxy 
formation. Larson (1974a,b) computed gaseous collapse models for ellipticals and noted 
that central star formation could continue until very late, fueled by dregs of gas falling into 
the center. He conjectured that central star density might depend on a delicate balance 
between the dissipation rate, global star formation efficiency, and mass loss via supernova- 
driven winds. 

Lauer (1983, 1985a) discussed the systematic properties of the centers of hot galaxies in 
light of new data. The high-density centers of small galaxies were consistent with Larson's 
gaseous infall picture, but the same theory predicted greater gas retention, and hence 
denser centers, in the deeper potential wells of bright galaxies, contrary to observations. 
An alternative scenario based on dissipationless merging also ran into difficulties because 
N-body simulations showed denser cores forming as galaxies successively merged (Farouki 
et al. 1983); this "departure from homology" has been confirmed with more modern N- 
body simulations (Barnes 1992; Makino & Ebisuzaki 1996). 
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Carlberg (1986, see also Ostriker 1980) used phase-space density arguments to con- 
strain the progenitors of hot galaxies. He noted that the high phase-space densities of small 
galaxies precluded their formation from purely stellar spiral disks because the phase-space 
density of spiral disks is low and phase-space density in dissipationless mergers remains 
constant or declines (Barnes 1992). However, HCM naturally incorporates gaseous merging 
as well as stellar merging (Schweizer 1986; Kormendy 1989; Kormendy & Sanders 1992), 
either during the main merger event (Negroponte & White 1983; Barnes & Hernquist 1991; 
Hernquist & Barnes 1991), during an earlier phase in which the disk progenitors develop 
bulges (Barnes 1992; Hernquist 1993), or during later gaseous infall. Thus it is no problem 
to create the high densities of small hot galaxies provided gaseous dissipation is present at 
some stage. 

The above papers considered mainly equal-mass mergers, but small satellites can 
also be accreted by dynamical friction (e.g., Tremaine 1976). Kormendy (1984) applied 
this concept to the capture of small, dense ellipticals by larger ones (see also Balcells 
& Quinn 1990) and predicted the photometric and kinematic signatures of such events: 
cores-within-cores, high surface-brightness centers, central velocity dispersion dips, central 
counter-rotation, and isophote twists. Many of these anomalies have since been found (de 
Zeeuw & Franx 1991; Barnes & Hernquist 1992), leading to the concept of (photometrically 
and/or kinematically) decoupled centers. Altogether, roughly a third of hot galaxies show 
such anomalies (de Zeeuw & Franx 1991), with core galaxies showing them roughly twice 
as often as power laws (Nieto et al. 1991b, Paper I). 

Hernquist & Barnes (1991) suggested an alternative way to make decoupled centers 
involving gaseous accretion and subsequent star formation in a cold inner disk. Decoupled 
subsystems are often dynamically colder than expected from purely stellar satellite accre- 
tion (Bender & Surma 1995, Franx & Illingworth 1988). Their stars are also stronger-lined 
than stars at the centers of small satellites, suggesting that local nucleosynthetic enrich- 
ment (and hence star formation) has taken place (Bender & Surma, 1988). 

To summarize, the consensus exists that gaseous dissipation plus in situ star formation 
are the key factors responsible for the high central densities of small hot galaxies. Both 
cores and power-laws probably also contain both stellar and gaseous material captured in 
accretions and mergers. Less clear is why central properties scale at all with galaxy mass 
and in particular why the centers of massive hot galaxies are so diffuse. 

7.2 Power Laws in Disky Galaxies 

We turn now to the relation between central and global properties. The basic question 
is whether the core/power-law distinction found for the inner parts is consistent with 
theories for forming the outer parts. 

The situation seems clearest for disky galaxies. The formation of disky hot galaxies, 
though not fully understood within HCM (Kormendy & Bender 1996), probably stems 
from the presence of significant quantities of gas during the latest merger(s). Several 
authors have noted that the high rotation, disky isophotes, and lack of minor-axis rotation 
of disky ellipticals imply global gaseous dissipation (Kormendy 1989; Nieto et al. 1991a; 
Bender et al. 1992). Recent simulations with gas (Barnes 1996, Barnes & Hernquist 1996) 



-17- 



show that just a few percent of the mass in gas is sufficient to destroy box orbits and 
impart high global rotation; the same gas can make disky isophotes if it forms stars. 

The steep central power laws of disky galaxies are plausible by-products of such 
gaseous mergers. Numerical models of gas-rich mergers (Barnes & Hernquist 1991, Hern- 
quist & Barnes 1991, Mihos & Hernquist 1994) have shown that angular- momentum trans- 
fer and dissipation can swiftly carry much of the gas in merging galaxies to the center 
of the remnant. Strong evidence for efficient gaseous infall is seen in nearby gas-rich 
merger remnants, which possess massive central clouds of gas (Sanders et al. 1988, Scoville 
et al. 1991) that may be fueling central starbursts (Bushouse 1987; Kennicutt et al. 1987; 
Leech et al. 1989; Condon et al. 1991; Solomon, Downes & Radford 1992; Kormendy & 
Sanders 1992). A major uncertainty is exactly where and how the stars form, and hence 
the shape and density of the resulting stellar profile; if anything, current models of star- 
burst cusps are too dense and compact (Mihos & Hernquist 1994). This problem aside, 
the high-density cusps of power-law galaxies seem broadly consistent with their formation 
in dissipative, gas-rich mergers 8 . 

7.3 Cores and Boxy Hot Galaxies 

The above discussion supports the notion that disky galaxies and their central power- 
laws were formed together in gas-rich mergers. Analogous arguments suggest that boxy 
galaxies formed in gas-poor mergers (Binney & Petrou 1985; Bender & Mollenhoff 1987; 
Nieto 1988; Nieto & Bender 1989; Nieto et al. 1991a; Bender et al. 1992). The distinctive 
shape, slow rotation, anisotropy, and minor-axis rotation of boxy galaxies are consistent 
with a large population of stars moving on box orbits in a triaxial potential created during 
a dissipationless merger (Barnes 1988, 1992). 

We therefore ask: are the core profiles of boxy galaxies simply the natural by-product 
of dissipationless stellar merging? To address this, the merging history of boxy galaxies 
can be simplified into two parts: an early phase in which centers originally formed, and 
a later phase involving the accretion of small satellite companions. Do core profiles form 
early, and do they survive later accretion? 

7.3.1 Early Formation 

The early formation of boxy, core galaxies is murky because their progenitors are 
poorly known 9 . Existing A-body simulations of dissipationless equal-mass mergers do not 

8 It may be significant that 6 out of 7 power-law galaxies in our sample with —20.5 > 
My > —22 are field objects, whereas 6 out of 9 core galaxies in the same range are in 
clusters. Kauffmann (1996), using Press-Schechter theory, has proposed that the clustering 
history of intermediate-magnitude ellipticals depends on environment — those in clusters 
are old, while many in the field are formed from recent mergers of gas-rich spirals. A 
preponderance of power-laws in field galaxies would be consistent with their formation in 
recent gas- rich mergers. In fact, two of the 6 field power-law galaxies are known merger 
remnants (NGC 596 and NGC 1700; Schweizer et al. 1990). 

9 It is clear, however, that bright ellipticals were not formed by simply merging to- 
day 's faint ellipticals. This is precluded by their very different stellar populations (Bender 
et al. 1992) — each type must have had its own progenitors. 



- 18 - 



develop cores — rather, pre-existing cores tend to shrink slightly due to non-homology, 
and the central density increases at each level of merging (Farouki et al. 1983; Barnes 
1992; Makino & Ebisuzaki 1996). Thus it appears that cores in luminous galaxies do 
not arise spontaneously in equal-mass merging, although the resolution of present N-body 
experiments is limited. 

The problem may be worse with wnequal-mass mergers, in which a smaller, denser 
component could sink to the middle, perpetuating a high-density center. This situation is 
discussed further under satellite accretion. On the other hand, progenitors of boxy galaxies 
may differ from today's hot galaxies and may not obey the same inverse correlation between 
mass and density. 

An entirely different way to generate low-density cores via stellar merging is to start 
with pure spiral disks. However, conventional density fluctuation spectra do not form 
spirals in the overdense environments that give rise to elliptical galaxies (Blumenthal 
et al. 1984, Bardeen et al. 1986). Yet a third way to make diffuse cores is via mass 
loss in stellar or AGN-driven winds, but the deeper potentials of luminous core galaxies 
should retain more gas, not less (Larson 1974b). A final possibility is to whip phase-space 
vacuum into centers during merging, for example via mergers of multiple subclumps (Weil 
& Hernquist 1996); however, it appears that the merging of the subclumps must be nearly 
simultaneous, which would be difficult to orchestrate for every core galaxy. 

7.3.2 Late Satellite Accretion 

The possible difficulty of forming cores may be matched or superseded by the even 
greater problem of maintaining them against satellite infall. In any merging hierarchy, the 
more luminous galaxies cannibalize the less luminous ones. It is plausible that the central 
region of the smaller galaxy will survive intact so long as its radius is smaller than the tidal 
radius imposed by the larger galaxy; this in turn implies that the regions of the smaller 
galaxy that are denser than the core of the large galaxy should survive. 

In fact, the centers of today's satellite galaxies are much denser than the centers of 
core galaxies (see Fig. 13, Paper I); typical power-law galaxies in the range My = —17.5 
to —22 are 100 to 1000 times denser at 10 pc than Abell brightest cluster galaxies (BCGs), 
and are 3 to 30 times denser at 100 pc, the inner boundary where accreted satellites can be 
detected by HST. The tidal argument therefore suggests that dense satellites should survive 
infall, filling in low-density cores of bright galaxies as proposed by Kormendy (1984). Yet 
every bright galaxy in our sample (except Fornax A) has a low-density core. 

To quantify this paradox, we estimate a typical satellite accretion rate for Abell BCGs. 
From counts of nearby companions and other data, Lauer (1988) deduced an accretion 
rate for BCGs of 0.2L* per Gyr, in close agreement with a theoretical estimate by Merritt 
(1985). There are 7 Abell BCGs in the current HST sample, all of which have large, low- 
density cores. 10 If satellite profiles are preserved during infall, accreted satellites over a 
particular magnitude range will be detectable. Small satellites have too little light, while 
large ones have profiles that are too similar to the BCG to make a difference. The profiles 

10 They are NGC 2832, NGC 4889, NGC 6166, NGC 7768, Abell 1020, Abell 1861, and 
Abell 2052. 
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in Fig. 1 imply 11 that satellites with My between —19.0 and —22.0 would be detectable in 
all 7 BCGs, and that those between —17.5 and —22.0 would be detectable in all but Abell 
2052. If BCGs have been accreting for 5 Gyr at the rate estimated by Lauer (1988), this 
translates to 2 detectable accreted satellites per BCG, or 13 total accretions in 7 galaxies. 

This estimate is conservative — gaseous accretion has been neglected, and the cur- 
rent accretion rate by BCGs is probably lower than average owing to the rise in cluster 
velocity dispersions with time. Restriction to Abell BCGs has excluded such near-BCGs 
as NGC 4874 in Coma (which was probably once the BCG of its subgroup), Virgo's BCG, 
and BCGs of smaller groups like Pegasus, Fornax, and Eridanus (where accretion is prob- 
ably faster owing to smaller velocity dispersions). Including such objects would double the 
number of primaries to 14 and raise the number of expected accretions to 26. However, 
no filled-in cores are seen in any of these BCGs. 

So far we have assumed that the inner portions of accreted satellites survive while 
sinking to the centers of their primaries, based on the tidal disruption argument. This ar- 
gument has been criticized by Weinberg (1994, 1997), who stresses that the time-dependent 
tidal force from the host galaxy can do work on resonant stars in a satellite galaxy even 
when the satellite is much denser than the host. Using semi-analytic perturbation theory 
and King-model profiles, Weinberg concludes that, if satellite and primary obey the global 
fundamental scaling law of Eq. (10), the satellite will be disrupted during its orbital decay 
if its mass exceeds 10 _3 -10 -2 of the primary mass. If this were true, low-density cores 
would remain unaffected by late accretion because of satellite disruption. 

The most relevant N-body simulation so far of satellite survival is a merger of two pure 
ellipticals with mass ratio 10:1 by Balcells & Quinn (1990). The density scaling of the small 
galaxy relative to the larger one approximately follows Eq. (10). At the end of the merger, 
the pre-existing core of the primary is filled in by an amount that would be detectable by 
EST. Further N-body models are in progress to check Weinberg's analytic results (Dubinski 
1997) and to simulate the dense central power-laws of real satellite galaxies (Minske & 
Richstone 1997). Realistic modeling of these dense centers may prove crucial. 

To summarize this section, the link between the centers of hot galaxies and their outer 
parts must be accounted for in the HCM picture. Many properties of both disky and boxy 
galaxies are naturally explained by appealing to a difference in the amount of gas present 
during the most recent merger(s). Disky galaxies, including their high central densities, 
suggest final mergers that were gas-rich. Analogous arguments concerning boxy galaxies 
are less clear: the global kinematics of these galaxies suggest final mergers that were gas- 
poor, but forming and preserving cores in such models may be difficult. An enlargement 
of the HCM model for core formation that includes BHs is considered in the next section. 

8. CORE CREATION BY MASSIVE CENTRAL BLACK HOLES 

11 Our criterion is that the net profile after infall be one magnitude brighter than 
presently observed at the inner resolution limit. This is sufficient either to erase a core in 
marginally resolved galaxies or create a tell-tale inner upturn in well-resolved cores. 



- 20 - 



High-resolution kinematic observations of galaxy centers strongly indicate that mas- 
sive BHs are normal constituents of the centers of hot galaxies (Kormendy & Richstone 
1995). Three new BH candidates have been discovered (Ferrarese et al. 1996 [NGC 4261]; 
Kormendy et al. 1997 [NGC 4486B]; Bower et al. 1997 [NGC 4374]), and the case for 5 
more has been strengthened (Harms et al. 1994 [M 87], Kormendy et al. 1996b [NGC 3115], 
Kormendy et al. 1996c [NGC 4594], Gebhardt et al. 1997 [NGC 3377], and van der Marel 
et al. 1997 [M 32]. BHs may play a key role in determining the central structure of galaxies, 
and no discussion of the central structure expected in HCM models would be complete 
without examining their influence. 

If both BHs and mergers are common among hot galaxies, two galaxies with pre- 
existing BHs will frequently merge. The BHs will spiral towards the center of the merger 
remnant, heating and perhaps ejecting the stars. This process may form the observed core 
in the merger remnant (Begelman et al. 1980; Ebisuzaki et al. 1991; Makino & Ebisuzaki 
1996; Quinlan 1997, Quinlan & Hernquist 1997). 

In what follows we assume that every hot galaxy contains a BH with average mass 

M. = 0.002 M gah (12) 

where M ga i is the mass of stars in the spheroid. The adopted coefficient 0.002 is a mean 
of estimates based on the energetics of AGNs and direct mass estimates of local BHs 
(see Appendix B). The assumption of proportionality in Eq. (12) is motivated by current 
data on local BHs (Kormendy & Richstone 1995; Kormendy et al. 1997), although the 
measurements show a scatter of at least an order of magnitude. 

The evolution of a pair of BHs in a merger remnant was first examined by Begelman 
et al. (1980). A recent comprehensive analysis is provided by Quinlan (1997). The two 
BHs are carried toward the center of the remnant by the general inward motion of the 
dense central parts during the merger but will not be exactly at the center at the end of 
the main merger phase. Subsequent migration of the BHs towards the center occurs on a 
slower timescale via dynamical friction from the background sea of stars. As the BH orbits 
decay, they form a bound binary BH whose semi-major axis a continues to shrink through 
dynamical friction. As the binary becomes more tightly bound, dynamical friction becomes 
less effective, and the characteristic decay time \d\ogaldt\~ 1 increases. Finally the binary 
orbit shrinks to the point that gravitational radiation or gas accretion dominates the decay, 
and rapid coalescence ensues. 

Decaying BHs lose most of their energy by heating the surrounding stars. The con- 
sequent puffing up of the galaxy was first examined by Ebisuzaki et al. (1991). Based 
on rough analytic arguments and iV-body models, they proposed that a merger of two 
galaxies with BHs would create a low-density core even if none previously existed. They 
argued that the mass of this core is approximately equal to the sum of the BH masses 
M. = mi + m2. If all galaxies start with the same ratio of BH mass to galaxy mass, 
this ratio would be unchanged by later merging, and thus the ratios rb/r e and M,/M ga i 
would remain constant. This was later seen in hierarchical merging N-body experiments 
with BHs (Makino & Ebisuzaki 1996) and also agrees approximately with observations 
(see below). 
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Quinlan and Hernquist have reexamined the evolution of binary BHs using scattering 
experiments and iV-body models (Quinlan 1997, Quinlan & Hernquist 1997). The following 
discussion is based on their results, which treat the infall of equal-mass pairs of BHs ranging 
in individual mass from 0.00125 to 0.04M ga i; the total BH mass M. = m\ + m-i ranges 
from 0.0025 to 0.08M ffa ;. Quinlan has kindly provided details of these models, which allow 
us to estimate the ratio of BH mass to core mass, a key quantity needed to compare to 
observations. 

The models start with a spherical galaxy whose density profile follows either a Hern- 
quist law, 

r (1 + r/r s ) 

or a modified Hernquist law with steeper slope in the inner parts: 

p( r ) °< riTTT / ^.5 - ( 14 ) 

r ±.o _|_ r /r s J 

The BHs are started on circular orbits at the half-mass radius r e of the galaxy. As the 
BH orbits shrink by dynamical friction, the stellar profiles develop cores. Final break radii 
r b were measured (by us) by locating the maximum of the logarithmic curvature of the 
projected mass surface density, which is equivalent to the definition in the nuker law. We 
define the indicative core mass as 

M core = Txr 2 b T,{r b ), (15) 

where E(rfo) is the projected surface density at rv The ratios M core /M., M core /M ga i, 
and rb/r e were tabulated for every model. M ga i is the stellar galaxy mass (dark matter is 
ignored) . 

The resultant core mass is approximately proportional to the BH mass but depends 
somewhat on the mean slope of the original mass profile over the region covered by the 
new core. For M,jM gal near 0.002 (Eq. 12) we find 

M core = (3.5 - 6.4) M„ (16) 

which translates to 

M core = (0.007 -0.012)M saJ (17) 

if M m /M ga i = 0.002. The range in parentheses reflects the two models in Eqs. (13) and 
(14). The scaling relation for break radii analogous to Eq. (17) is found to be 

r h = (0.02 - 0.06)r e . (18) 

These results imply that the orbital decay of a BH creates an indicative core mass 
that is 3-6 times the BH mass. This is larger than the core mass M core ~ M. estimated by 
Ebisuzaki et al. (1991), and larger than the ejected mass M e j £ 2M. found by Quinlan & 
Hernquist. Our explanation is that the ejection of a given mass can create the impression 
of a more massive core simply due to the precise definition of r core as the point of maximum 
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logarithmic curvature in the profile — the exact location of r core depends sensitively on 
how it is defined. The larger indicative core mass in Eq. (16) comes closer to matching 
the observed indicative core mass for M 87, which is ~ 10M,. 12 

The theoretical predictions of Eqs. (17) and (18) are compared to observed core 
luminosities and break radii in Figs. 9a and 9b. Indicative core luminosity for observed 
galaxies is defined (analogously to M core ) as L core = Tvr^Ib and is computed from the core 
parameters in Table 2. The dashed lines are power-law fits to the observed data derived 
by assuming unit log slope and weighting all points equally. The observed relations are 

L core = 0.012 L gah (19) 

and 

r h = 0.03 r e . (20) 

The gray areas represent the ranges covered by the theoretical predictions in Eqs. (17) 
and (18); in plotting Fig. 9a, it is assumed that M core /M ga i = L core /L ga i, which should 
be true provided (M/L)v for the stars does not vary strongly between 10" and r e . 

Observed core radii are within the range predicted by the models for M./M ga [ = 0.002, 
while core luminosities are near the upper boundary of the predicted range. The predicted 
trends as a function of luminosity are generally matched, although observed core masses 
may increase as a steeper-than-unity power of total mass (cf. Eq. 8). 

The Quinlan-Hernquist models confirm the suggestion (Ebisuzaki et al. 1991, Makino 
& Ebisuzaki 1996) that mergers of galaxies containing massive BHs can generate cores 
with roughly the size and luminosity indicated by the observations. The presence of BHs 
in cores also helps to defend cores against accretion: small dense satellites accreted by a 
primary core galaxy may be disrupted by the BH before they sink to the center. 

Despite these encouraging results, models of core formation by massive BHs remain 
uncertain. In particular, (1) the existing simulations do not yet explore the full range 
of BH mass ratios and initial conditions appropriate for merging galaxies; and (2) initial 
BH formation has simply been posited in ad hoc fashion in all hot galaxy progenitors (see 
Haehnelt & Rees 1993). 

There are also further problems to be considered: 

(i) Why do core profiles exhibit weak cusps? Perhaps the slow shrinkage of the BH binary 
naturally forms a cusp, either because the stars are flung into elongated orbits by the 
binary, or in the same way that cusps are formed when a single central BH grows 
adiabatically (Peebles 1972; Young 1980; Quinlan, Hernquist & Sigurdsson 1995). 
However, the N-body models of Quinlan & Hernquist do not show such cusps at 
present resolution. Alternatively, gas infall into a core previously formed by a BH 
binary might steepen the profile to create a cusp (Begelman et al. 1980; Young 1980). 
Gas is apparently collecting now at the centers of at least some core galaxies (e.g., 
M87, Ford et al. 1994). 

12 Based on the core parameters of M 87 in Table 2, the fitted nuker-law global M/Ly 
value for M 87 of 10.2, and M. = 3 xlO 9 M (Harms et al. 1994, as scaled by Kormendy 
& Richstone 1995). 
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(ii) What is the relation of nuclei to BHs? Do they signal BHs, compete with BHs, or 
possibly feed BHs? The nucleus in NGC 3115 (Kormendy et al. 1996b) has a stellar 
mass of ~ 3 x 10 7 M crammed into a tiny volume of radius ~2 pc around a BH that 
is 50 times more massive. The stellar density approaches 1O 6 M pc -3 , and typical 
orbital velocities exceed 1000 km s -1 . It is a puzzle how stars could have formed 
in such an environment, where gas clouds are likely to be colliding at high velocities 
while bathed by intense radiation from the BH. Perhaps the nucleus and BH formed 
in different progenitors that later merged. 

(hi) Our discussion so far has stressed correlations of core and global properties. In fact 
there are outliers such as Fornax A, which has a very small core for its luminosity 
(Fig. 4a). Fornax A is peculiar and is probably still in the throes of a major merger 
(Schweizer 1980). Perhaps the inner regions have not yet settled down to their final 
state, giving us a clue to the time scales involved in core scouring; or gas (as signaled 
by the copious dust) may be (re) forming a stellar cusp, although there is no sign of 
young stars in the color map (Shaya et al. 1996). 

(iv) BHs appear to be associated with the hot component of spiral galaxies; late-type 
spirals such as M33 have little or no central BH (Kormendy & McClure 1993). When 
two late-type galaxies merge, they are believed to form an elliptical, but this will not 
have the central BH that is required for core formation in subsequent hierarchical 
merging. 

(v) Most important, if cores are formed by merging binary BHs, why do power-law galaxies 
at intermediate magnitudes (—20.5 > My > —22) not have cores with size as given 
by Eq. (18)? BHs appear to be just as common in power-law galaxies (Kormendy & 
Richstone 1995). Perhaps power laws can be regenerated by star formation from fresh 
gas supplied by the latest merger. However, to avoid being ejected by the BH binary, 
the new stars must form after the BH binary shrinks, which poses a timing problem 
if BHs sink to the center more slowly than gas. 

9. SUMMARY AND CONCLUSIONS 

We have assembled inner surface-brightness profiles for 61 dynamically hot galaxies 
available in the HST archive as of June 1993. Fits of the nuker law (Eq. 2) to deconvolved 
profiles are used to compute values of break radius r& and break surface brightness 1^. 
These are supplemented with ground-based data from the literature on rotation, isophote 
shape, and velocity dispersion. These data are used to produce updated versions of the 
central parameter diagrams for hot galaxies. 

The inner surface-brightness profiles of hot galaxies can be divided into two types as 
discussed in Paper I. Core galaxies have a sharp knee or bend in the profile, akin to the 
analytic cores of King models or the isothermal sphere but with a shallow cusp at small 
radii. Power-law galaxies have profiles that are steep and rather featureless in log-log 
coordinates with no detectable core at . 1 resolution. 
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Cores appear only in galaxies brighter than My = —20.5 (Fig. 4a); core size and 
luminosity are roughly proportional to galaxy size and luminosity (Eqs. 18 and 19). Power- 
law galaxies are fainter than My = —22. In the overlap region from My = —20.5 to —22, 
the two types coexist and profile morphologies vary widely — upper limits to core size in 
some power-law galaxies are at least 100 times smaller than the core sizes of other galaxies 
at the same luminosity. 

The scatter in central properties in the overlap region correlates with global structure: 
core galaxies tend to be boxy and rotate slowly, while power-law galaxies are disky and 
rotate rapidly. Preliminary evidence suggests a further correlation with environment in 
that core galaxies tend to be found in dense groups and clusters while bright power laws 
are preferentially found in the field. 

Cores populate a fundamental plane (FP) that is analogous to and roughly parallel 
to the global FP for elliptical galaxies. The scatter about this plane (in logn,) is 0.12 
dex, about 30% larger than the analogous scatter about the global FP. Some of this extra 
scatter may come from massive BHs, which may inflate central velocity dispersions in some 
galaxies. 

A set of self-consistent scaling relations for core galaxies is presented that expresses 
core size, density, and other quantities as a function of Ly (Eqs. 3-10). These scaling laws 
are projections of the FP. A major conclusion is that small hot galaxies are much denser 
than large ones, by a factor of up to 1000 at a radius of 10 pc. 

The last part of the paper attempts to relate the central parameter relations of hot 
galaxies to the process of galaxy formation and evolution. We suggest that the presence of 
dense power-law centers, disky isophotes, and rapid rotation in low-luminosity galaxies all 
point to their formation via dissipative, gas-rich mergers. The analogous arguments about 
core galaxies are less clear: the boxy isophote shape and slow rotation of these luminous 
objects suggest formation by dissipationless mergers, but cores may be difficult to form 
and maintain in such events. For example, core galaxies seem at present to be accreting 
small dense satellites in sufficient numbers to fill in their low-density cores, at least if the 
satellites survive their orbital decay to the center, an issue that is still in dispute. 

We explore an alternative model for core formation based on merging BHs. The 
model assumes that BH binaries are formed in galaxy mergers; the binary orbit decays 
by dynamical friction, ejecting stars from the center of the merger remnant, enlarging any 
previous core, and scouring out a new one where none existed. Simulations of this process 
by Quinlan & Hernquist (1997) yield a reasonable match to the radii and masses of observed 
cores if every hot galaxy contains a central BH of average mass M. = 0.002 M ga \. This 
value for M, is consistent with BH mass estimates in AGNs and local BHs (see Appendix 
B). Whether or not BHs are the dominant agent in creating cores, their role in shaping 
the central structure of hot galaxies is likely to be significant if they are as common and 
as massive as recent estimates suggest. 

The main goal of this paper is to explore systematic trends in the central structure 
of hot galaxies and the possible relations between their present central structure and their 
formation history. By strengthening the link between the central structures of hot galaxies 
and their global properties such as luminosity, shape, and rotation, HST has helped to 
open an important new window on galaxy formation. 
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APPENDIX A: THE NATURE OF POWER LAWS AND THE 
FREQUENCY OF INNER DISKS IN POWER-LAW GALAXIES 



Jaffe et al. (1994) and Ferrarese et al. (1994), like us, divide hot galaxies into two types 
based on inner surface-brightness profile. Those called by us cores with a strong break and 
low central surface brightness they term Type I, and those called by us power-laws with 
no break and high central surface they call Type II. There is no discrepancy between us 
as to division into classes based on profile shape. 

Jaffe et al. go on to identify core galaxies in a general way with slowly rotating boxy 
galaxies, and power-law galaxies with rotating disky galaxies. This distinction resembles 
ours but differs in important details. For example, Jaffe et al. envision that the centers of 
hot galaxies either have or do not have small inner disks. Such disks seen edge-on are what 
create the high surface brightness power-law profiles of Type II galaxies. These inner disks 
are furthermore associated with the global, outer disks of rotating disky galaxies. It is 
the frequent association between inner disks and outer disks that creates the link between 
power-law profiles and disky rotating galaxies in their picture. 

Jaffe et al. believe that the high surface brightness profiles of power-law galaxies are 
produced only when an inner disk is seen edge-on. Specifically they state: "Most of the 
characteristics. . .that discriminate Type I from Type II are explained by disk components 
seen at high inclination angles. For example, the higher central surface brightness in Type 
II systems is caused by the nuclear disk seen close to edge-on. . .If one of these [disky] 
galaxies [i.e., a Type II] were viewed face on, it would appear much more like a Type I 
galaxy." Thus, their view is that observed profile type is due to a combination of intrinsic 
properties plus viewing aspect. All power-laws have inner disks - they are the disks that 
happen to be seen edge-on. Core-type profiles on the other hand are a mixture; many are 
intrinsic cores that lack inner disks, while some fraction are disks seen face-on that are 
masquerading as cores. 

This interpretation of cores versus power laws differs importantly from our own. Our 
view is that the difference between core and power-law profiles is intrinsic to the hot stellar 
component and has no direct connection with a disk, whether seen edge-on or face-on. 
Power laws remain power laws at any viewing angle, as do cores. Paper I presented initial 
arguments against the edge-on disk explanation for power laws. We have since undertaken 
a more comprehensive comparison with power-law galaxies in the present data set. Briefly, 
we find that the small sample of Virgo galaxies analyzed by Jaffe et al. was abnormally 
dominated by late- type edge-on SO galaxies. Nearly every power-law galaxy they detected 
was such an object. A high frequency of edge-on inner disks in such a sample is therefore 
understandable. Our sample is larger and contains many power-law ellipticals that are not 
flattened and show no sign of either an inner or an outer disk. This and other evidence to 
be described leads us to conclude that power-laws are independent of inner disks and are 
thus a feature of the hot stellar component alone. 

We carefully examined deconvolved V-band images of all 61 galaxies in the present 
sample. Thirteen of the 14 Jaffe et al. galaxies were available in the archive, and we 
looked at all of them. Seven of these were admitted into our sample (they are included in 
Table 1). The remaining six galaxies were rejected for the following reasons: (1) too much 
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dust to derive a reliable surface-brightness profile or class the object as a core or power 
law (NGC 4261, NGC 4342, NGC 4374, NGC 4476); (2) a potential double nucleus and 
unclassifiable profile (NGC 4473); (3) no clear spheroidal component (NGC 4550); and (4) 
interfering spiral arms (NGC 4476). Compared to our sample, the Jaffe et al. power-law 
galaxies are much later in type, diskier, more edge-on, and more subject to dust and other 
peculiarities that potentially interfere with reliable measurement of the spheroid profile. 

If the edge-on disk interpretation were correct, then all or most power-law galaxies in 
our sample should show evidence of disks. To test this, we inspected each image within 
a 10" radius for an inner disk comprised of either stars or dust. We also looked for spiral 
arms, which we took as another indicator of a disk. We visually assessed isophote shape 
as a function of radius, checking for changes and comparing to measured values of a^. We 
tried to correlate changes in ellipticity and shape with kinks or "ledges" in the brightness 
profile — such correlations might signify the edge of a disk. We also subtracted the profile 
fits given in Paper I and looked for signs of a residual disk. Table 4 summarizes the results 
of this visual inspection. On the basis of this evidence, we assigned a final score to each 
galaxy indicating the likelihood of an inner disk. The values are (no sign of a disk), 1 
(possible disk), 2 (probable disk), and 3 (disk plainly visible). These scores are given in 
the table, along with comments. 

The results of these efforts are summarized in Fig. 10, which plots inner disk score 
versus ellipticity for power-law galaxies. The hypothesis that all power-law objects have 
edge-on inner disks seems unlikely because: (1) Half the sample shows little or no evidence 
of any inner disk, including several highly flattened objects. (2) Power-law galaxies are 
not concentrated at high ellipticities, in contrast to the sample of Jaffe et al. (3) Most 
objects with inner disks are bulges that also have outer disks (types SO, Sa, or Sb), a 
point also made by Jaffe et al. However, this weakens the case for ubiquitous inner disks 
in power-law ellipticals because they lack such outer disks. (4) One power-law galaxy, 
NGC 3599, shows face-on spiral structure, showing conclusively that it cannot be edge-on. 
(5) If known bulges are excluded, the ellipticity distribution of the remaining power-law 
ellipticals is nearly the same as that of core ellipticals (which lie within the rectangle in 
Fig. 10 but are not plotted individually), and neither type shows much evidence for inner 
disks. Thus the evidence for inner disks is weak in both power-law and core ellipticals. 

A less restrictive hypothesis, not put forward by Jaffe et al., is that power laws are 
associated with a high surface-brightness inner disk, period, whether seen either edge-on or 
face-on. This also seems unlikely because there are several flattened power-law ellipticals 
that must be close to edge-on yet show no sign of a disk (point 1 above). The kinematic 
properties of M 31 and M 32 are also relevant here. Both of these would be typical power 
laws if seen at a distance, yet both are hot and slowly rotating at radii of a few arcsec, the 
claimed size of inner disks in other power-law galaxies. In neither galaxy is there any hint 
from kinematics that the high central surface brightness is associated with a disk. 

To summarize, present evidence does not favor the ubiquitous presence of high surface- 
brightness inner disks in power-law galaxies, whether edge-on or not, though such disks 
are certainly present in some cases. Rather we believe that the difference between core 
and power-law profiles more probably reflects an intrinsic difference in the spheroidal light 
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distribution between the two types. Finally we note cautiously that one of our low-disk- 
score galaxies, NGC 3377 (disk score = 1), has since revealed a dust disk in recent WFPC2 
images. Final conclusions about the frequency of inner disks in ellipticals and bulges should 
therefore await a new body of high-quality WFPC2 images. 



APPENDIX B: THE MEAN BH MASS PER HOT GALAXY 

The following argument adapted from Tremaine (1997) summarizes the evidence for 
the frequency and masses of BHs in the centers of hot galaxies. 

The integrated comoving energy density in quasar light (as emitted) is (Chokshi & 
Turner 1992) 

u = 1.3 x 10" 15 erg cm" 3 , (21) 

independent of H and Q. If this energy is produced by burning fuel with an assumed 
efficiency e = AE/ (AMc 2 ), then the mean mass density of dead quasars must be at least 
(Soltan 1982, Chokshi & Turner 1992) 



u _ _ K /0.1 
ec 



- - 2.2 x 10 5 ( — )M Mpc" 3 , (22) 



assuming that the Universe is homogeneous and transparent. 
The mass of a dead quasar may be written 



0.1 



M. = ^= 7Xl 0X( I ^)( I ^)^j, (23) 

where Lq is the quasar luminosity and r is its lifetime. An upper limit to the lifetime 
is the evolution timescale for the quasar population as a whole, ~ 10 9 y; however, upper 
limits to BH masses in nearby galaxies and direct estimates of the BH masses in AGNs 
both suggest that the actual masses and lifetimes are smaller by a factor 10-100 (Haehnelt 
& Rees 1993), which implies M, = 10 7 -10 8 M Q . 

To focus discussion, we adopt a strawman model in which a fraction / of all galaxies 
contains a central BH, and BH mass is proportional to galaxy luminosity. Thus M. = XL, 
where T is the (black hole) to (galaxy) mass-to-light ratio. The luminosity density of 
galaxies is j = 1.5 x 10 8 L Q Mpc -3 in the blue band (Efstathiou et al. 1988, adjusted to 
our Hubble constant of H = 80 km s _1 Mpc _1 ). Thus 

This value is an average over the light of all local galaxies. However, if we assume that 
massive BHs are found chiefly in the centers of hot galaxies (Kormendy & Richstone 1995), 
the above number can be converted to the BH mass-to-light ratio per hot component by 
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noting that approximately 30% of the local B-band light is emitted by such components 
(Schechter & Dressier 1987). Correcting for this and converting to the V band yields 



*-T(t)5-. 

where T v is now the estimated ratio M,/Ly per hot component and fh is the fraction of 
hot galaxies with BHs. 

A second estimate of T v from quasars comes from dividing the typical dead quasar 
mass derived above, M. 1O 7 5 M , by the typical luminosity of a bright hot component, 
8.5 x 1O 9 L (Binggeli, Sandage & Tammann 1988, adjusted to H = 80 km s^Mpc" 1 ). 
The result is T y ss 0.004. Consistency with Eq. (25) then requires fh ~ 1 if e « 0.1, or 
that most or all hot galaxies contain a BH. 

A final method for estimating T v uses individual BH masses for 6 moderately well 
established BHs in nearby hot galaxies 13 (Kormendy & Richstone 1995, Table 1). Using 
an estimate of global stellar M/ L v based on nuker-law fits as described in Table 3 14 yields 
the logarithmic mean value T v = 0.016. This is in reasonable agreement with T v = 0.004 
from quasars in view of the likelihood that these best BH candidates are more massive 
than average. 

For further discussion, we assume that every hot galaxy contains a BH and adopt for 
T v the logarithmic mean of the quasar and BH values: 

T v = 0.008. (26) 

The corresponding value of M,/M ga i is then 

M./M gal = 0.002, (27) 

where M ga \ comes from Ly(M/L)y and (M/L)y is the above-mentioned mean global 
mass-to-light ratio for the 6 candidate BH galaxies. This is the value of M,/M ga i per hot 
component adopted in Section 8. 



13 The Galaxy and NGC 4258 are omitted for lack of accurate luminosities of their hot 
components, and NGC 3115 uses the new value of M. = 2 x 10 9 M Q from Kormendy 
et al. (1996b). 

14 For the present purpose, the mass fits were renormalized to fit a at 10" using the ratio 
R a = cr (10") /ao from Table 2. This was done to avoid possible contamination of <jq by a 
BH. The logarithmically averaged (M/Ly) for the 6 galaxies is 4.0 in solar units (Hq = 80 
km s _1 Mpc _1 ). 
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FIGURES 



Figure 1: V-band surface-brightness profiles of 55 ellipticals and bulges from HST. All 
were observed in the WFPC1 Planetary Camera through filter F555W and were decon- 
volved using the Lucy-Richardson algorithm as described in Paper I. Core galaxies (see 
Section 2) are plotted as solid lines, and power-law galaxies are plotted as dashed lines. 
"Mean radius" is the geometric mean of the semi-major and semi-minor axes of the isopho- 
tal ellipse. 

Figure 2: HST surface-brightness profiles of M 31 and M 32, as seen locally and near 
Virgo (24 times farther). To simulate Virgo, the nearby profile was binned by a factor of 
24, convolved with the WFPC1 point-spread function, and deconvolved with 80 iterations 
of the Lucy-Richardson algorithm. 

Figure 3: Division of the sample into cores and power laws. The figure plots logarithmic 
inner slope of the surface-brightness profile, 7, versus angular break radius, 9^, from fits to 
the nuker law (Eq. 2). Galaxies with log 9^ > —0.8 are well resolved and divide into two 
groups with high and low 7. Dashed lines connecting the near and far versions of M 31 
and M 32 indicate potential resolution effects on other power-law galaxies. A galaxy must 
have 7 < 0.3 and a well-resolved break radius to be classed as a core. Galaxies within the 
box comprise the "Core" sample. 

Figure 4: HST measurements of central parameters of hot galaxies, as a function of 
absolute V magnitude. Hubble type and nucleus types are taken from Table 1; "bulges" 
are SO-Sb galaxies, and nt for power laws are limits rf m and \i\ im from Table 2. M 31 
and M 32 are plotted twice: asterisks show data as observed, and tails indicate their 
positions as they would appear 24 times farther away near Virgo. The small black square 
is the SO galaxy NGC 524, which is the only core within a bulge. The apparent turndown 
in surface brightness at faint magnitudes in panel (c) is probably a resolution effect (cf. 
M 32). Effective radii are plotted in panel (d), to be compared with break radii in panel 
(a): the strong impressions of scatter at intermediate magnitudes (—22 < My < —20.5) 
and of two types of galaxies in panel (a) are absent in panel (d). 

Figure 5: Break radius r& versus distance. The dashed line is the adopted dividing line 
for cores in Fig. 3 (log 9b = —0.8). Above this line, a core-type profile will be seen as a 
resolved core, below it will be classed as a power-law. The trend versus distance is opposite 
to what one would have expected if cores and power-laws were merely an artifact of angular 
resolution — core galaxies are on average more distant than power-laws. Moreover, most 
of the sample is close to Virgo in distance (log-D ~ 3.2), yet contains both cores and power 
laws, confirming that the two types are intrinsically different. 

Figure 6: Various densities at radius .1 plotted against absolute magnitude. Mass 
densities are derived by normalizing nuker-law surface- brightness fits to central oq. The 
symbols are the same as in Fig. 4. Model details are given in the text and notes to Table 
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3. Panel (a) luminosity density; panel (b) mass density; panel (c) peak Maxwellian phase- 
space density. Note the range of almost 10 6 in density in all three panels. Turndowns for 
small galaxies are probably an artifact of resolution (cf. M 32). 

Figure 7a: ReplotofFig. 4a with symbols indicating rotation speed (i>/cr)*. Slow rotators 
(filled symbols) have (v/cr)* < 0.51; fast rotators (open circles) have (v/a)* > 0.51. Bulges 
lacking data are classed as fast rotators. Galaxies with core profiles are indicated by the 
enclosing squares; all others are power laws. The data indicate a tendency for fast rotators 
to have power-law profiles. 

Figure 7b: Same as Fig. 7a but with symbols indicating isophotal shape a 4 /a. Galaxies 
are classed as disky if a 4 /a > 0.4, otherwise as boxy/neutral. Irregular profiles with 
variable a 4 /a are also classed as boxy/neutral. Bulges (Hubble types SO-Sb) are classed 
as disky. The data indicate a tendency for disky galaxies to have power-law profiles. 

Figure 8: HST measurements of central parameters of hot galaxies in fundamental-plane 
space. Symbols are the same as in Fig. 4. Tails on M 31 and M 32 (asterisks) show the 
effect of moving these galaxies 24 times further away to the vicinity of Virgo. Resolution 
effects on other power-law galaxies may be similar and are indicated schematically by the 
limit flags. Panel (d) shows the fundamental plane rotated about the do axis and viewed 
edge-on (for cores). The rotation chosen uses the same power-law combination of r& and 
h used for the global fundamental plane of elliptical galaxies by Faber et al. (1987) and 
is consistent with their best core plane within the statistical errors. The rms scatter about 
the central plane (cores only) is 0.12 dex, which is 50% greater than the scatter about the 
global plane. This increase may be due in part to the influence of central BHs on uq. Five 
BH candidates from Kormendy & Richstone (1995) are marked in panel (d). 

Figure 9: Core versus global properties. Panel (a) plots indicative core magnitude (com- 
puted from L core = Kr\lb) versus total magnitude. The dashed line is a mean fit assuming 
unit logarithmic slope (see text). Panel (b) is similar but compares break radius r& to 
effective radius r e . The shaded areas represent predictions of decaying BH binary models 
(Quinlan & Hernquist 1997). 

Figure 10: Inner disk prominence versus ellipticity for power-law galaxies. Disk score 
(Table 4) is a visual estimate of the evidence for an inner disk: = no evidence, 1 = slight, 
2 = probable, 3 = definite. Symbols are the same as for power laws in Fig. 4. Core galaxies 
are not plotted individually; they lie within the rectangle. If bulges (SO-Sb) are ignored, 
there is little remaining tendency for power-law galaxies to be highly flattened, as might 
be expected if they were due to inner disks seen edge-on (Jaffe et al. 1994). There is also 
little tendency for inner disks to appear in flattened galaxies, which would be expected if 
they were aligned with the outer isophotes. See Appendix B for further discussion. 
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Fig. 1 : Inner Profiles of 55 Ellipticals and Bulges 
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Fig. 2: M31 and M32 as Seen Locally and at Virgo 
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Table 3. Local Quantities from Nuker-law Fits. 
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Notes to Table 3. 

Col. 1: Name (see Table 1). 

Col. 2: Hubble type (see Table 1). 

Col. 3: Absolute V mag of bulge component from Table 1. 

Col. 4: Profile class: fl = core; \ = power law. 

Columns 5-11 are based on nuker law fits to the brightness profile from Table 2 (see text). M/L is assumed constant, and the velocity dispersion 
is isotropic everywhere. The models are mass-normalized by fitting to the central dispersion, (To, assuming that it is measured through a 2"x 2" 
aperture (in perfect seeing). If (To is missing, columns 5, 7, 9, 10, 11, 15, and 16 are blank. 

Col. 5: Normalized M / Ly in solar units. 

Col. 6: j(0 .1) is the luminosity density at a radius of . 1 in Lq pc — 3 . 

Col. 7: Equivalent Maxwellian peak phase-space density at .1, given by / o/(27T(T 2 ) 3 / 2 , in units of Mq pc 3 (km s 1 ) 3 . 

Col. 8: V-band luminosity within a sphere of radius .1, in Lq. 

Col. 9: Mass within a sphere of radius .1, in Mq. 

Col. 10: Projected line-of-sight velocity dispersion at .1, in km s — 1 . 

Col. 11: Local relaxation time in yr at a radius of . 1 from Equation (8-71) of Binney and Tremaine (1987). The Coulomb log has been set to 
log(0.4M(< O'.'l)/O.7M ), where M(< O'.'i) is from column 9. 

Col. 12: Physical radius corresponding to .1, in pc. 

Col. 13: Luminosity density as in column 5, but at a radius of 10 pc. 

Col. 14: Enclosed luminosity as in column 8, but within a sphere of 10 pc. 

Col. 15: Radial velocity dispersion at 10 pc (note: not line-of-sight as in column 10). 

Col. 16: Local relaxation time as in column 11, but at 10 pc. 



Figure 3: Division into Cores and Power-laws 
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Table 4. Inner disk survey of power-law galaxies. 
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SO galaxy. Subtracted model shows thin disk close to edge-on. 
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Inner disk score: = no sign of a disk; 1 = possible disk; 2 = probable disk; and 3 = disk plainly visible. 



Fig. 4. Central Parameters vs. Absolute Magnitude 
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Fig. 5: Cores and Power-Laws vs. Distance 
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Fig. 6. Densities at 0.1" vs. Absolute Magnitude 
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Fig. 7. Correlation with Rotation and Isophote Shape 
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Fig. 8. Projections of the Core Fundamental Plane 
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Fig. 9. Core vs. Global Properties 
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Fig. 10: Inner Disk Prominence vs. Ellipticity for Power-Laws Galaxies 
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Notes to Table 1 . 



Col. 1: Name. "A" objects are first-brightest Abell cluster galaxies. "VCC" objects are dwarf Virgo 
Cluster E's from Binggeli, Sandage, and Tammann (1985). M 31 = NGC 224; M 32 = NGC 221; 
M 87 = NGC 4486. The entries for NGC 221 and NGC 224 followed by "V" denote M 31 and 
M 32 as seen 24 times further away, just beyond the Virgo Cluster. 

Col. 2: Hubble type from the Second Reference Catalog of Bright Galaxies (de Vaucouleurs et al. 1976, 
RC2): 1=E, 2=E/S0 or SO, 3=Sa or Sb, 4=dE or dSph. 

Col. 3: Group number from Faber et al. (1989). 

Col. 4: Adopted distance in units of km s — 1 . Group membership information is used where avail- 
able. The following distance estimates are combined in descending order of weight: 1) Surface- 
brightness fluctuation method (Tonry, unpublished). 2) D n — a average for Faber et al. (1989) 
groups with 3 or more members. 3) Radial velocity (group or single) corrected by smooth pecu- 
liar velocity field as determined by POTENT (Dekel et al. , in preparation). 4) Radial velocity 
in cosmic microwave background frame (for distant objects). 

Col. 5: Distance in column 4 converted to Mpc using H = 80 km s -1 Mpc" 1 . 

Col. 6: Fully corrected B9^ of bulge component, from Faber et al. (1989) where available (E's only, 
magqual >2), otherwise mostly from the RC2 and/or the Third Reference Catalog of Bright 
Galaxies (de Vaucouleurs et al. 1991, RC3). Special sources of magnitudes: VCC ellipticals: 
Binggeli and Cameron (1993); Bender (priv. comm.). Abell brightest cluster galaxies (BCGs): 
add 0.21 mag to VI mag of Hoessel et al. (1980), where the constant is determined from 7 
galaxies with VI and B9^ in common. 

Bulge/total brightness ratios (in B-band) for SO's and spirals: NCG 224 = 0.24 from de Vau- 
couleurs (1958). NGC 524 = 0.88 and NGC 1023 = 0.51 from Kormendy and Illingworth 
(1983). NGC 2814 = 0.24 from Boroson (1981). NCG 3115 = 0.94 from Capaccioli et al. (1987). 
NGC 3384 = 0.59 from Burstein (1979). NGC 4594 = 0.93 from Burkhead (1986). NGC 7332 
= 0.58 from average of SO's in Simien and de Vaucouleurs (1986). NGC 7457 = 0.38 from 
Burstein (1979) and Kormendy (1977). 

Col. 7: B-band Galactic extinction Ag from Faber et al. (1989) or RC3. 

Col. 8: Fully corrected (B — V) from Faber et al. (1989) where available, otherwise from RC3. Bulge 
(B — V) for some SO-Sb's is a guess. 

Col. 9: Profile class: fl = core; \ = power law. 

Col. 10: Degree of nucleation: ++ = severe, + = moderate. M 31-in- Virgo is the dividing line between 
severe and moderate (see Figure 2). 

Col. 11: Upper limit to break 9^ (in arcsec) for power-law galaxies only. Based on simulated model 
profiles convolved with HST PSF and then deconvolved to match data. Most values are from 
Paper I, but some are new here. Used in preference to 9^ for power-law galaxies, see text. 

Col. 12: Lower limit to break surface brightness /j,^, for power law galaxies only. Goes with 0^ m in 

column 11. Raw value uncorrected for Galactic extinction in V mag arcsec -2 . Most values are 
from Paper I, but some are new here. Used in preference to ^ for power-law galaxies, see text. 

Col. 13: Break radius 9^ in arcsec from nuker-law fit. Most fits are from Paper II, but some are new 
here. Used for core galaxies, see text. 

Col. 14: Break surface brightness /i^ in V mag arcsec -2 from nuker-law fit. Raw value uncorrected for 
Galactic extinction, in V mag arcsec -2 . Used for core galaxies, see text. 

Col. 15-17: a, [3, and 7 from nuker-law fit, mostly from Paper II. 

Col. 18: Sources of HST images: 1 = Paper I; 2 = WFPC1 Team GTO program (Ajhar et al., in 
preparation); 3 = WFPC1 Team GTO program, unpublished; 4 = Lauer et al. (1992a); 5 = 
Lauer et al. (1993); 6 = Lauer et al. (1992b); 7 = Lauer et al. (1991); 8 = Shaya et al. (1996); 
9 = Grillmair et al. (1994); 10 = Jaffe et al. (1994; Virgo survey); 11 = Forbes et al. (1995; 
kinematically decoupled cores). 



Notes to Table 2. 



Col. 1: Name (see Table 1). 

Col. 2: Hubble type (see Table 1). 

Col. 3: Absolute V mag of bulge component based on Bj^, (B — V) , and distance from Table 1 (Hq 
= 80 km s — 1 Mpc -1 ). No K-correction or cosmological corrections have been applied. 

Col. 4-5: Logarithm of central velocity dispersion, <jq, and ratio of central dispersion to projected dis- 
persion at 10", <to/<tio. The geometric mean of a\Q on major and minor axes is used if both 
are available. <jq is assumed to equal to the rms luminosity-weighted value through a 2" x 2" 
aperture in l"FWHM seeing. Quoted values are weighted means from the following sources: 
Bender & Nieto (1990); Bender et al. (1994); Bertola et al. (1988); Binney et al. (1990); Davies 
& Birkinshaw (1988); Davies & Illingworth (1983); Dressier & Richstone (1988); Dressier & 
Richstone (1990); Efstathiou et al. (1980); Efstathiou et al. (1982); Faber et al. (1989); Fisher 
et al. (1995); Franx et al. (1989a); Fried & Illingworth (1994); Gonzalez (1993); Jedrzejewski & 
Schechter (1988); Jedrzejewski & Schechter (1989); Kormendy (1982b); Kormendy (1988); Ko- 
rmendy & Illingworth (1983); Kormendy & Richstone (1992); Kormendy & Westpfahl (1989); 
Scorza & Bender (1995); Tonry (1984); van der Marel et al. (1994); Whitmore et al. (1985); 
Young et al. (1978). 

Col. 6: Dimensionless ratio (i>/(t)* = (v/a)/(v/a) n a t e (see Davies et al. (1983) for definition). Most 
values have been taken from the literature; those calculated by us assume (v / <j) bi a te = — 
e)]i fromBender (1988). Published sources of (v/a)*: Bender et al. (1992); Bender et al. (1994); 
Bosma et al. (1985); Davies et al. (1983); Fisher et al. (1995); Gonzalez (1993); Jedrzejewski & 
Schechter (1989); van der Marel (1991). 

Col. 7-8: Ellipticity and isophotal shape parameter, a^/a (xlOO), as defined by Bender & Mollenhoff 
(1987). a^/a is the coefficient of the cosine distortion term expressed as a percentage of major 
axis length. C4 is also the cosine term but expressed as a percentage of surface brightness 
fluctuation around an elliptical isophote. a^/a is related to C4 by the local brightness gradient 
as follows: 



a^/a = C4 



dlogl 



d\ogr 



-1 



(b/a 



,1/2 



where the second term corrects from the unit circle to the major axis (Bender et al. 1988). 
Where this formula is needed, we have simply taken | (/log// dlogr | = 2 and ignored the second 
term. In estimating a^/a from profiles, we have followed Bender et al. (1988) and used an 
average between 10" and 60". In the text, galaxies with a^/a > 0.4 are classed as "disky" while 
all others (including irregular 04 /a 's) are classed as "boxy". 

Sources for e and a^/a: Bender et al. (1989); Faber et al. (1989); Franx et al. (1989b); Goudfrooij 
et al. (1994); Jarvis & Freeman (1985); Kent (1983); Lauer (1985b); Lugger et al. (1992); Nieto 
et al. (1991c); Peletier et al. (1990); RC2; van den Bosch et al. (1994). For disk galaxies, we 
have attempted to find e for the bulge only, but those for NGC 1023, NGC 2841, NGC 3384, 
and NGC 7332 are global axial ratios from the RC2. 

Col. 9: Profile class: fl = core; \ = power law. 

Col. 10/12: Logarithm of break radii, in pc. Angular values from Table 1 (6^ tm and 6^) have been converted 
to linear values (r^ m and r^) using the distance in Mpc from Table 1. No cosmological curvature 
corrections have been applied. 

Col. 11/13: Break surface brightnesses ^ m and /j,q from Table 1 corrected for Galactic extinction (but not 
for K-correction or cosmological dimming) . 

Col. 14-16: a, [3, and 7 from nuker-law fits, repeated from Table 1. 

Col. 17-18: Logarithm of the effective radius, in pc, and effective V-band surface brightness from Faber 
et al. (1989; NGC 1700 from Gonzalez 1993). Hq = 80 km s" 1 Mpc" 1 is assumed, and no 
cosmological curvature corrections have been applied. \i e is the mean V-band surface brightness 
within r e and is calculated from the B-band value in Faber et al. (1989) using (B — V) from 
Table 1. 



